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ABSTRACT 



We study the fluctuations of standard thin accretion disks by linear analysis 
of the time-dependent energy equation together with the vertical hydrostatic 
equilibriu m and the equation of state. We show that some of the simulation 
results in lHirose et al.l (J2009bJ), such as the time delay, the relationship of power 
spectra, and the correlation between magnetic energy and radiation energy, can 
be well understood by our analytic results. 



Subject headings: accretion, accretion disks - hydrodynamics - instabilities - MHD 



Introduction 



Aperiodic X-ray fluctuations have been observed from both galactic black hole binaries 
(BHBs) and active galactic nuclei (AGNs) (Uttley et al. 2005; McHardy et al. 2006). The 
Power Spectral Density (PSD) of such variability is generally modeled with a power law, 
P(f) oc /~*, where P(f) is the power at frequency /, and varies with frequency. In 
the soft state, the PSDs of both BHBs and AGNs have a steep slope with ~ 2 at high 
frequencies, flatting to a shallow slo pe with ~ 1 b elow a bend frequency /b, which is 
typically around 10Hz for BHBs (see iKing et al.l 12004 and references therein). The PSDs 



in the hard state are more complex. The origin of variability is not well understood yet. 
However, it is highly tempting to relate this variability to the magnetohydrodynamic (MHD) 
turbulence, which is believed to drive the accretion process (IBalbus fc Hawleylll99ll ). Some 
works followed this path through numerical simulat ion, and typically used proxies for the 
radia tion rather than a direct measure of luminosity (IHawley fc Krolikl|200ll ; iNoble fc Krolik 
20091 ). It remains uncertain whether the proxies for radiation are appropriate to describe the 
luminosity fluctuations. In the present work, we will show that the variability of magnetic 
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energy (stress) of the standard thin disk is different from that of radiation for short time-scale 
fluctuations. 



Recently, shearing box simulations of stratified magnetorotational turbulence ( IHirose et al. 



2009bj) showed that fluctuations in the magnetic energy (stress) lead those in the radiation 
energy with roughly a thermal time-scale, and a correlation is found between the stress and 
total pressure. Moreover, the disk is found to be thermally stable, which is, howe ver, in 
conflict with the disk theory ( iLightman fc Eardleyl Il974j ; IShakura fc Sunyaevl 1 1 9 73h . The 
discrepancy reveals that the correlation found in the simulation may be different from the 
a-prescription. For example, such a correlation may be related to the energy equation or 
result from the feed back from pressure to stress, which is not in the form of the standard 
a-prescription (e.g.. lLin et al.ll2011t ICiesielski et al.ll2012l ). Since the dissipation of magnetic 
energy will heat the gas of accretion flow, the perturbations in the magnetic energy will pro- 
duce corresponding fluctuations in the internal energy and therefore in the pressure. Then, 
there should exist a correlation and delay between the stress and the pressure. In the present 
work, based on the energy equation, we will investigate the relationship between fluctuations 
of the viscous heating and the inducing fluctuations of the radiative cooling. 

The paper is organized as follows. The relationship of fluctuations of the viscous heating 
and the inducing fluctuations of the radiative cooling is derived in Section 2. A comparison 
of analysis and simulation is presented in Section 3. Conclusions and discussion are made in 
Section 4. 



2. Radiative cooling fluctuations induced by viscous heating fluctuations 



2.1. Energy equation 



In the context of standard thin accretion disk (IShakura fc Sunyaevill973l ). t he vertically 
integ rated energy equation in cylindrical coordinates (r, <f), z) takes the form (e.g-. lKato et al. 
20081 ): 

dE ,„ _ dln£ „«91n# + _ 



+ n- 



dt dt 

where H is the vertical height of the disk, and £ (= 2pH) and II(= 2pH) are the surface 
density and the vertically integrated pressure, respectively. The gas internal energy E and 
the radiative cooling rate Q^d P er un ^ area are expressed as 



E = R 



rad 



E 



gas 



3(1-/3) + 



g 

7-1 



n 



(2) 



-3 - 



16acT 4 

Qrad = ~W ' (3) 
where E Tax ± and E gas are respectively the internal energy of radiation and gas, (3 is defined as 
the ratio of the gas to the total pressure, i.e., (3 = n gas /n, 7 is the ratio of specific heating, and 
T is the temperature on the equatorial plane of the disk. The opacity R is generally dominated 
by the electron scattering (k 6S ) in radiation-pressure-dominated accretion disks, where R can 
be regarded as a constant. On the other hand, if the opacity is dominated by the free-free 
absorption (Kg), R will vary with the temperature and the density. The viscous heating Q^ is is 
due to the dissipation of magnetic energy and turbulent kinetic energy in magn etoturbulent 



disks , and is dominated by the dissipation of magnetic energy in simulations (j Simon et al. 



2009|). 



For the fluctuations with a time-scale less than the viscous time-scale, the variation of 
£ can be neglected, and Equation ([1]) is therefore simplified as 

dE BlnH ^_ 

In order to study the induced fluctuations of radiative cooling, we adopt the vertical hydro- 
static equilibrium 

tf K H> = 5, (5) 
and the equation of state, which can be approximately expressed as 

n = n gas + n rad = -^st + \qT^e. (6) 



2.2. Relationship of fluctuations 

We use the subscripts "0" and "1" to describe the unperturbed and perturbed quantities, 
respectively. We would stress that, the amplitude of fluctuat ions is assumed to b e small in 
our linear analysis. In simulations (e.g., Figures 3 and 4 of iHirose et al.ll2009bl ). however, 
the amplitudes can be significantly large. Nevertheless, the linear analysis may reveal the 
relationship of fluctuations of physical quantities. Combining Equations ©-© with R = n es , 
we have + 

QvisA 



Q 



vis,0 



A 7 ! 
T ' 



(7) 



where the dimensionless parameter A is expressed as 
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(4 -3/3) ( 7 -l) 
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and the thermal time-scale £th takes the from (with Q 



vis,0 



3an fi K /2): 



Eq 



'th 
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rad,0 
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vis,0 
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2/3 



3( 7 



(9) 



We choose 7 = 5/3 and a = 0.02 (e.g., iHirose et al.ll2009al ) for numerical calculations. The 
variation of A with /3 is shown by the solid line in Figure [TJ The other two parameters, 
A s (dashed line) and A rad (dotted line), will be introduced by Equations f[T3"j) and ffTEj) . 
respectively. 

By assuming that the time-dependent component of fluctuations takes the form of 
exp(iut), e.g., Q^ is 1 /Q^ is oc exp(icot), we have the following relationship from Equation ([7]): 



To 



Q 

iAut th + 1\Q 



vis,l 



+ 

vis,0 



(10) 



where (Qtis^/Qtisfi)^ an d (^Ti/Tq)^ represent the fluctuations with uj of the viscous heating 
and those of the radiative cooling, respectively. We would stress that Equation f lTOj) is a key 
relationship in the present work. 



3. Comparison of analysis and simulation 
3.1. Time delay between magnetic energy and radiation energy 



Equation ( ITOj) can be modified as 



Q 



vis,l 



yj (Autth, 



2 _|_ l \ ^?vis,0 



exp {—iut 



del J 



TT 



where the delay time i<jei of the radiative cooling compared with the viscous heating takes 
the form: 



tdei = — arctan (Aut 



th) 



(12) 



Obviously, this equation implies tdei ~ ^th for long time-scale fluctuations. 



Figure M shows the variatio n of tdei f° r three different values of /3. The simulations 
for (3 ~ 0.1 ( jHirose et al.ll2009bl ) showed that fluctuations of magnetic energy lead those of 
radiation energy by 5 — 15 orbit periods (t rb), roughly a thermal time. Moreover, Figure 5 of 
Hirose et al.l (l2009bl ) indicates that significant variability occurs in the range 0.01 < /-t or b ^ 
0.1, which is equivalent to 0.01 < u/Qk < 0.1, corresponding to the region between the two 
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vertical dot-dashed lines in Figure 2. As shown by the solid line, the delay in our analysis is 
around 3 — 17t or b, which is consistent with the simulations. In addition, we would point out 
that the delay between viscous heating and magnetic energy 
is negligible compared with the thermal time-scale. 



0.5t orb , iHirose et al.ll2009bf ) 



Furthermore, simulations have been done for the gas-pressure-d ominated case (jHirose et al 



20061 ) and the case that gas and radiation pressures are comparable (IKrolik et al.ll2007l ). The 
delay in those simulations is ~ 2t orb for ~ 0.8 and ~ 5t or b for ~ 0.5. In our analysis, 
as shown by the dotted and dashed lines, the delay is around 2 — 3t or b for = 0.8 and 
2 — 12t or b for = 0.5 in the range 0.01 < oj/Qk ^ 0.1, which is again co nsistent with the 
simulations. Note that R in the simulations for ~ 0.8 (IHirose et al.l 120061 ) is dominated by 
the free-free absorption. In such case, for a simple approach, we modify the parameter A as 
Aff by considering the free-free absorption instead of the electron scattering in Equation (3): 



(4 -3/3) ( 7 -l) 



(11.5 + 4.5/3) [0 + 3(1-/3) (7 



7-60 + 



20 7(3 (4 - fry) (1 - 0) 
7-1 (7 - 1) (4 - 30) 



(13) 



The profile of As is shown by the dashed line in Figure 1. 



3.2. Power spectrum relationship 

In this subsection, we will show a comparison between our analytic normalized power 
spectrum of radiation energy P r ^ d and that in simulations P ra d- Equation (ITT]) provides 
the relationship of the power spectrum between the radiative cooling P coo i and the viscous 
heating P vis : 

P cooi(f) = : ; 7z * ^ vis(/)- (14) 



1 + (27rAft thi 



2 ' 

A 



Then the analytic power spectrum of volume-integrated radiation energy P^ d is expressed 
as 2 

^rtd(/) = — 2 P vis(/), (15) 

radV ' 1 + (27iAft th ) 2 K J V ' 

where the quantity A ra ^ is derived from Equations (J2J) , (]5j) , and (EJ: 

. = (Prad.l/Prad.o) _ - 4 - 30 

A rad _ (4Ti/To) i + 4 (1 + ^- ^> 

The variation of A rSi d with is shown by the dotted line in Figure [TJ Since the viscous 
heating is mainly due to the dissipation of magnetic energy, it is plausible to have 

P vis (/) « \(V2nft dis ) 2 + ll P B (f), (17) 
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where tdis is the dissipation time-scale of magnetic energy, and Pb is the power spectrum of 
magnetic energy. The explanation for this relationship is presented in Appendix A. Then, 
Equations (fT5|) and (|T7|) provide an analytic relationship between the power spectrum of 
radiation energy and that of magnetic energy: 



pA 



rad 



(f) 



A 2 



(V2nft dis ) 



1 + (2tt Aft 



-PbU). 



tbJ 



The simulations (Equation (17) and Figure 5 of iHirose et al.ll2009bl ) showed the profiles 
of P B and P rad : 

8.7 x 10- 6 / -113 , 

10-7-3-65, 

6.1 x 10" 9 /- 2 - 38 , 
2.3 x 10- 10 /- 3 - 91 , 



PbU) 



Prsdif) 



f< 0.171, 

/> 0.171, 

/ < 0.118, 
/ > 0.118. 



(19) 



(20) 



In our Figure 3, we replot the above Pb and P ra d with the dotted and dashed lines, respec- 
tively. In addition, according to Equation (fT8"|) . we plot the analytic power spectrum P 1 



rati 



with the solid line. The values of Pb(/) i n Equation ( ITBl are taken from the above simula- 
tion results (the dotted line). The parameters for calculating P^, are A = 2.1, v4 ra d = 1.8 
(corresponding to (3 = 0.1), and tdis = 0.5i or b ( IHirose et al.ll2009bl ). As shown by the solid 
and dashed lines, our analytic P r ^ d agrees well with P ra d in simulations. 



3.3. Correlation between magnetic energy and radiation energy 

The energy equation implies a correlation between the viscous heating (magnetic energy) 
and the pressure. Based on Equation (j7]), the correlation can be read as: 

S^4^. (21) 

Owning to the same reason, there should be a correlation between viscous heating and 
magnetic energy, i.e., Qti S ,i/ Qtis,o — Eb,i/Eb,o- With Equation (16), the above equation can 
be modified as 

Eba „ ^E^^ (22) 

Eb,0 ^-rad Prad,0 

As shown in Figure 1, there exists A rad = 1.8 for (3 = 0.1, so we have the relationship: 

Eb oc P r ° ad 5 , (23) 



which is close to the correlation found in simulations, e.g., E B oc P 9 ^ 1 (IHirose et al.ll2009bl ). 
The difference in the index may be related to the following two reasons: (1) the analysis 
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is quite simple, particularly in dealing with the vertical radiative cooling; (2) the feedback 
from pressure to stress makes significant contribution. 



It is worthy to note that iHirose et al.l (l2009bl ) also provided an explanation for the 
correlation with a toy model based on the energy equation. By using the correlation between 
£th an d E TSu d (tth oc E° ad ) obtained in simulations, they derived the correlation between Eb 
and _E ra d (Eb oc E t ^). In our analysis, we choose Q~ d to replace their radiative cooling 
term E^/t^. We will show below that our results are quite similar to theirs. 



Equations (j3J) and ffTB"]) can provide the relationship: 

Q- rx E^ AmA 

Vrad ^ "^rad 



In simulations (IHirose et al.ll2009bl ). the thermal time is calculated by 



t 



th 



E/Q 



rad' 



With Equations (J2J) , fl5J), (jbj) and (j24j) . the above equation can be reduced to 

£ th oc £ r s ad , 



(24) 



(25) 



(26) 



where 



3(1 - /3)(4 - 3/3)( 7 - 1) - 3/3(1 + /3) 



(8 + [/3 + 3( 7 -l) (1-0)] 
For radiation-pressure-dominated accretion flows, Equation (|2"7|) can be simplified as s' 



(27) 



1/A-ad- Thus, Equation ( 1221) indi cates the relationsh ip £?b oc E^ , which is consistent 



with the correlation in the toy model (IHirose et al.ll2009bl ). Moreover, in the case of (3 = 0.1, 
Equation ( 127|) gives s 1 = 0.41, thus t t h ^rad 1 ) which is close to the correlation found in 
simulations, e.g., t t h oc E®^? in simulation 1112a and t t h oc E®^ in simulation 1126b of 



Hirose et al.l (l2009bh . 



4. Conclusions and Discussion 

In the present work, we have studied the fluctuations of standard thin disks by linear 
analysis of the time- dependent energy equation together with the vertical hydrostatic equi- 
librium and the equation of state. Our analytic results show that the delay between magnetic 
energy and radiation energy is consistent with that in previous simulations. In addition, the 
analytic power spectrum of radiation energy agrees well with that in simulations. Moreover, 
the correlation between magnetic energy and radiation energy can be well understood by the 
analysis, with an index (0.55) being close to that in simulations (0.71). 
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As indicated by Equation (j!4p . there may exist a break frequency /b r ~ l/(2irAt t h) in 
Pcooi- The frequency /br m &y be associated wi th the high-freq uency break observed in the 
power spectra of luminosity fluctuations (e.g., iMcHardyl 120101 ) . since its value in the inner 
region of disk is close to that of observed high frequency break. Moreover, the difference 
between P coo i and P v j S for / > /br, shown by Equation (TH|) . shoul d be taken into ac count 



in modeling th e high-frequency variability of quasar luminosity (e.g., iMushotzky et al.l 12011 



Zu et al.ll2012l ). In addition, the similar frequency break may also occur in radiatively i neffi- 



cient accre tion flows, such as advect ion-dominated accretion flows (INarayan fc Yilll994i ) and 
slim disks (lAbramowicz et al.lll988l ). 
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A. The relationship between P vis and P B 

In this Appendix, we try to derive the relationship between P v - m (f) and P B {f) as shown 
by Equation (JTTj). The evolution of magnetic energy E B {t) can be simply described as 



dE B (t) 
dt 



G B {t)-D B (t), 



(Al) 



where G B (t) and D B (t) are respectively the generation and dissipation rate of magnetic 
energy. With small amplitude perturbations in Equation (lAlj) . we have 



iujt 



dis 



E B ,i 
E B n 



Gb,i 
Gb,o 
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B.l 



D 



B,0 



(A2) 



or 



(ut dis ) 2 P B (uj) = P g {uj)+P d {^ 1 Gin 
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B,0 
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B.l 
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B,0 
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BA 
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B,0 



(A3) 



where D B0 = G B:0 , t^is = E B0 /D Bi0 , the symbol "*" represents the complex conjugate 
number, the power spectrum of G B and D B are 
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B,l 



G 



Bfi 



Pn(uj) 



D 



BA 



Dbm 



-9 - 



Magnetic fields in the acc retion disk present expo nential growth owing to the magne- 
torotational instability (MRI, iBalbus fc Hawleyl Il99l[ ). followed by the dissipation due to 
some destructive mechanisms. Since the rise and decay phases of channel modes are similar 
(e.g., Figure 4 of ISimon et al.l 120091 ) . it is plausible to believe that P G and P D is compara- 
ble for oj < 7r/tdis- With the consideration of the delay (~ t$s) between Gb and Db, the 
relationship between Gb and Db can be modeled as 



D 



B,l 



G 



B,l 



exp (— iujtdh 



<Db,o J u \Gb,o, 
Substituting this relationship into Equation ( 1A3[) . we obtain 



(A4) 



Pn(oo) 



{cot, 



dis 



Pr(uj) 



2 (1 — coswtdis) 
The above relationship is applicable for oj < n/t^is- 



(cot 



dis 



Pb(J). 



(A5) 



The fluctuations with oo > tt /t<n s in and those in Db is unclear, and thus it remains 
uncertain for Pg(co) and Pd(oj). However, it may be plausible to believe that Gb and Db 
are decoupled with each other for uj > 7r/t dis . If we further assume that Pg(co) ~ Pd(cu), 
Equation flA3j) can be reduced to 



Pd(co) 



(utdh 



Pb(co). 



(A6) 



We use this equation to describe the relationship of Pd{co) and Pb(co) for oj > n/tdis- Based 
on Equations (1A5I) and (1A6j) . a general form of relationship between Pjj and Pb may be 
simply described as 

P D (f) « \(V27rft dis ) 2 + ll P B (f). (A7) 



Since the turbulent kinetic energy follows the fluctuating magnetic energy ( iHirose et al. 
2009bh and magnetic dissipation dominates over kinetic dissipation, we obtain 



Pvisif) « PdU) ~ \(V2nft dis ) 2 + ll PbU) 



(A8) 



which is the exact form of Equation ( TT7|) . It should be noted that the relationship between 
P vis and Pb in the short time-scale range (co > 7r/t dis ) is tentatively used in the present work. 
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Fig. I. — Variations of A (solid line), A Tad (dotted line), and A s (dashed line) with (3. 
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Fig. 2. — The analytic delay of radiation energy compared with magnetic energy for /3 — 0.1 
(solid line), 0.5 (dashed line), and 0.8 (dotted line). The vertical dot-dashed lines represent 
two specific frequencies cu/f2 K = 0.01 and 0.1. 
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Fig. 3. — A comparison of P r ^ d (solid line) and P ra d (dashed line), where P r ^ d is calculated 
with Pb (dotted line). 



